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^ ■ A relativistic single particle model is used to calculate the inclusive (e, e') reaction 
O ■ 

. from A =12, 40, 56, 197, and 208 nuclei in the quasielastic region. We have shown 

^ ' that this model provides a very good description of the available experimental cross 

Q"^ ' sections when they are dominated by the quasielastic process. In this paper we 

i"^ ' use this model to investigate the dependence of y-scaling on electron kinematics, 

I , particularly the electron scattering angle, for a range of squared four momentum 

O ■ 

^ I transfer 0.20 — 0.80 (GeV/c) . In this kinematic domain. Coulomb distortion of 

the electron does not significantly affect scaling, but final state interactions of the 
, knocked out nucleon do affect scaling particularly when the nucleons have lower 

•r:^ energies. In general, we find that scaling works for this reaction, but at lower values 

o ■ 

' of the four momentum transfer, the scaling function does have some dependence on 

■ the electron scattering angle. We also consider a modification of y-scaling to include 

, small binding energy effects as a function of Z and A and show that there is some 

^ ' improvement in scaling. 

H 

; PACS numbers: 25.30. Fj 

Medium and high energy electron scattering has long been acknowledged as a useful tool 
to study nuclear structure and nuclear properties, especially in the quasielastic region where 
the process of knocking out nucleons is dominant. Many experiments 

mm 

have been 

performed on medium and heavy nuclei at incident electron energies less than 1 GeV where 
astic processes can be avoided. There are also a number of theoretical 



contributions of ine 

works 0,y,0,i8i,y,Q,[ii, 



12 



As we have noted before 



which have been compared to the measured responses. 
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12l |. the Fermi gas model in the impulse approximation 
roughly describes the inclusive (e, e') cross sections, but fails to provide a good description 
of the structure functions. A good theoretical description of quasielastic scattering requires 
two ingredients before one can compare experimental (e, e') data from light to heavy nuclei 



to theory. One of them is a model for the nuclear transition current and the other is some 
provision for the inclusion of electron Coulomb distortion effects for medium and heavy 
nuclei. In the early 1990's, the Ohio University group 5l,[l^,[l^ treated the electron Coulomb 
distortion exactly for the reactions (e, e') and (e, e'p) in the quasielastic region using partial 
wave expansions of the electron wave functions in the distorted wave Born approximation 
(DWBA). The Madrid group subsequently reported a similar analysis for the exclusive 
(e, e'p) reaction. However, the DWBA calculations do not allow a separation of the cross 
section into a longitudinal part and a transverse part and are numerically challenging, and 



computational time increases rapidly with 
avoid these difficulties, Kim and Wright [10 



ligher incident electron energies. In order to 
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IGj developed an approximate treatment of 



the electron Coulomb distortion which does allow the separation of the cross section into a 
longitudinal part and a transverse part. 

In addition, we found a model which provides a very good description of quasielastic 
scattering processes from nuclei for both the inclusive and exclusive cases. It is a relativistic 
single particle model which requires the wave functions of bound and continuum nucleons 
and a transition current operator. The bound state wave functions are obtained from solving 
the Dirac equation in the presence of the strong vector and scalar potentials [17] and for the 
inclusive (e, e') reaction where the knocked out nucleons are not observed, the continuum 
wave functions are solutions to a real potential so as not to lose any flux. At low energies, this 
potential is just the same as the bound nucleon potential and thereby guarantees current 
conservation and gauge invariance. However, it is known from elastic proton scattering 
that the continuum potential becomes weaker with increasing proton energy. Therefore, for 
higher energy processes we use a nucleon potential whose strength has been fltted to proton 
elastic scattering 18|, ll9|. In a recent paper 20|], we found excellent agreement between our 
model with the higher energy experimental data from SLAC [21] for the quasielastic (e, e') 
scattering on ^^C, ^^Fe, and ^^^Au at the squared four momentum transfer of approximately 
0.20 - 0.30 (GeV/c)^ by using the energy- dependent real potentials (which are weaker than 
the bound state potentials) for the outgoing nucleons. Note that these calculations do not 
conserve nucleon transition current, but we calculate all four components of the transition 
current in order to minimize errors due to lack of current conservation. Note that the Madrid 



group 



151 ] has used a very similar relativistic model for (e, e'p) reactions. The excellent 



agreement of this model with experimental data allow us to use it as a tool for investigating 



y-scaling of the quasielastic components of the inclusive (e, e') even in kinematic regions 
where inelastic contributions have large contributions. 

n n n 

Since the pioneering work by West there have been many experimental [4, |23| and 



theoretical 
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29j investigations of y-scaling from nuclei. In this paper 



we propose to use our relativistic mean field single particle model with the inclusion of 
Coulomb corrections to investigate the approach to y-scaling at intermediate values of the 
four momentum transfer. Using ?/-scaling, the measured cross section for the inclusive (e, e') 
reaction can be written as a product of the electron-nucleon cross section times a function F 
which is related to the momentum distribution of nucleons in the nucleus and is a function 
of momentum transfer q and energy transfer uj. For the case of large momentum transfer, 
ihe function F should depend only on a single variable y which is a function of uj and q 
261]. Scaling is expected to be valid for the very large momentum transfer region, but it 
may be broken by final state interactions in the quasielastic region and/or electron Coulomb 



distortion effects. A SLAC experiment 



2l| was performed at the squared four momentum 



transfer of 0.23 - 2.52 (GeV/c)^ and the data at y < exhibited scaling at large Q^. 
Recently, there have been additional experimental data from JLAB js^ at 4.045 GeV, 
but scaling of the quasielastic process cannot be demonstrated since the energy transfer 
is sufficiently high that pion production is a significant contribution to the cross section. 
Note that both of these experiments have been carried out at relatively forward electron 
scattering angles. Furthermore, we cannot compare our model predictions to these data 
since we have not yet included inelastic processes (meson production, etc.) in our model. 
However, as noted above we can use our model to investigate the scaling of the quasielastic 
contributions to the cross section. 

Using non-relativistic models, the authors in Ref. I27I analyzed y-scaling of the quasielas- 
tic electron scattering in few-body system, complex nuclei, and nuclear matter. Within the 
framework of the plane wave impulse approximation (PWIA), they investigated the effects 
of the final state interaction, the binding correction, and the nucleon-nucleon correlations. 
They pointed out that the relation between the scaling function and the momentum dis- 
tribution does not exist at finite momentum transfer because of the final state interaction 
and the binding correction. More recently, the Madrid group in collaboration with Donnelly 



29|, 



31 



32| have investigated scaling using a semi-relativistic model and note that the strong 



scalar and vector potentials in the final state of relativistic models result in a breakdown of 



scaling and result in different scaling functions for longitudinal and transverse responses. 

In this paper, we initially calculate ^/-scaling at the squared four momentum of ap- 
)roximately 0.2 - 0.3 (GeV/c)^ comparing with the experimental data measured at SLAG 



2l| for ^^C, ^^Fe, and ^^''Au, and Bates |3|] for ^°Ca, where the quasielastic contribution is 
kinematically isolated from pion production. We then investigate the effects of final state 
interactions in our relativistic model and electron Coulomb distortion on scaling. Since we 
are considering cases with large outgoing nucleon energies, we do use an energy dependent 
final state interaction. Finally, we introduce a new y-scaling variable in order to solve the 
non-scaling problem in the presence of the final state interaction from different target nuclei 
in the same four momentum transfer range. 

In the plane wave Born approximation (PWBA), where the electron wave functions are 
described by the Dirac plane waves, the cross section for the inclusive quasielastic (e, e') 
reaction is written as 



dujdVtp 



where q^^ = tu^ — = —Q^ is the squared four momentum transfer, aM = 
(a/2i?)^[cos^(6'e/2)/sin^(6'e/2)] is the Mott cross section, and the longitudinal and trans- 
verse structure functions which depend on the three momentum transfer q and the energy 
transfer u are Sl and St- Explicitly, the structure functions for a given bound state with 
angular momentum jf, are given by 

with the outgoing nucleon density of states pp = The z-axis is taken to be along the 

momentum transfer q and the z-components of the angular momentum of the bound and 
continuum state nucleons are fib and Sp, respectively. The Fourier transform of the nuclear 
current J^(r) is simply given by 

m = J J^(r)e''i"ciV, (4) 

where J'^(r) denotes the nucleon transition current. The continuity equation could be used to 
eliminate the 2;-component (A^^) via the equation A^^ = — ^A'o if the current is conserved, but 



+ (tan^ I - ^)5T(g,^)} , (1) 



since we use an energy dependent final state interaction [20| we violate current conservation 
and to minimize errors calculate A^^ directly. The nucleon transition current in the relativistic 
single particle model is given by 

J^{t) = e^p{v)J^M^) , (5) 

where is a free nucleon current operator, and ipp and ipb are the wave functions of the 
knocked out nucleon and the bound state, respectively. For a free nucleon, the operator 
comprises the Dirac contribution and the contribution of an anomalous magnetic moment 
/i-r given by J'^ = Fi(gp7^ + F2{qf^)j^a'^''q^. The form factors Fi and F2 are related to the 

2 

electric and magnetic Sachs form factors given by G^; = Fi + ^^-^2 and Gm = Fi + ^tF2 
which are assumed to take the following standard form: 

= ^2 — = 7 7-TY 5 (6) 

(1 - %Y i^^T + 1) 

where the standard value for is 0.71 (GeV/c)^. 

The y-scaling function is defined as the ratio of the measured cross section to the off-shell 
electron-nucleon cross section as follows: 

where a^p {(Jen) denotes the off-shell electron-proton (neutron) cross section acci defined by 



Ref. 



33|. In Eq. ([7]), Z and are the number of protons and neutrons, and M is the mass 



of nucleon. The scaling variable y 3J] is given by 



uj + MA = {M^ + q^ + y^ + 2yqf'^ + {Ma-i^ + y^Y'\ (8) 

where Ma is the mass of the target nucleus and Ma-i is the mass of the ground state of the 
A — 1 nucleus. The point ?/ = corresponds approximately to the peak of the quasielastic 
scattering and y < (y > 0) corresponds to the small (large) u region. 

In Fig. [H we calculate the y-scaling functions for ^^C, '*°Ca, ^^Fe, ^^^Au, and ^osp]^ 
neglecting both the final state interaction of the exiting nucleons and the electron Coulomb 
distortion. This is normally referred to as the plane wave impulse approximation (PWIA). 
The squared four momentum transfer is approximately 0.2 (GeV/c)^ < < 0.3 (GeV/c)^. 
The scaling function for all of these function of the scaling variable y are very 

similar-with deviations from the mean less than 10% which is probably due to binding 



energy effects. Note from tlie caption that all of these examples were calculated with electron 
scattering angles of 45° since, as we will show later, there is some breakdown of scaling at 
lower energies when the electron scattering angle is changed significantly. 

In Fig. [21, we show the same results as in Fig. [H but with electron Coulomb distortion 
turned on. We refer to this case as the distorted wave impulse approximation (DWIA). For 
heavier nuclei, electron Coulomb distortion shifts the scaling curves toward the right side. 
However, the deviations from scaling due to Coulomb distortions do not seem to be too 
large and the largest effect occurs for larger positive values of y. In Fig. [3l we show the 
same results as in Fig. [H but with both the final state interaction and the electron Coulomb 
distortion included. 

Clearly the final state interaction in this kinematic range leads to rather large violations 
of scaling which is also observed experimentally. For example. Fig. H] shows the comparison 
of our theoretical results with the experimental data for ^^C, ^^Fe, and ^^''Au measured 
at SLAC 2l|. The incident electron energy is E = 2.02 GeV and the electron scattering 



angle is ^ = 15". Under these kinematic conditions, the energy transfer is below the pion 
production peak (except for large positive values of ?/)and hence the quasielastic peak is 
well separated from inelastic processes. Note that our model, see also Q, describes the 
quasielastic process for all three nuclei quite well and furthermore, scaling is observed for 
the two heavier nuclei (^^Fe and ^^^Au) when the quasielastic process dominates {y < 0.1). 
However, the lighter nucleus (^^C ) does not scale with the heavier ones. This deviation is 
due to a combination of electron Coulomb distortion and final state interaction effects. 

In order to improve scaling we introduce a modified scaling variable y' given by y' = 
y — {N / Z)\Eb\av where \Eb\av is the average of the absolute binding energy for all the bound 
nucleons. This approach is meant to remove differences among nuclei including binding 



energy effects as suggested in Ref. [26|, |27| and may permit scaling for quasielastic scattering 



from different nuclei. In Fig. [5], we show the new scaling functions for ^''^Pb, ^^^Au, ^^Fe, 
^°Ca, and ^^C at the same kinematics as Fig. [31 Comparing Figs. [3] and [5] clearly shows 
that y' improves scaling to some degree although some deviations from scaling remain. 

In Fig. [HI we compare the scaling function to the experimental data measured at Bates 
[sl on '^°Ca for three cases. The solid line and • are for incident electron energy E = 739 
MeV and scattering angle 9 = 45.5°, the dotted line and ■ for E = 372 MeV and 9 = 90°, 
and the dashed line and ▲ represent E = 367 MeV and 9 = 140°. The scaling functions 



have the same shape and the peaks he at the same position, but the backward angle results 
(dashed curve and A) have a different magnitude from the others. From these calculations, 
we see that scaling at these somewhat lower energies is affected by the electron scattering 



angle. This result is in agreement with the studies reported by Caballero et al. [32] that 
final state interactions affect longitudinal and transverse contributions differentially since 
changing the electron scattering angle changes the longitudinal and transverse contributions 
to the cross section. Note that these calculations agree with the experimental data relatively 
well as in our previous results , |35|. In order to investigate the dependence of scaling at 



lower four momentum transfer values on 6, we calculated y-scaling with Coulomb distortion 
and final state interactions for only the protons in ^°^Pb (Z=82 and N=0) as shown in Fig. 
[7] and for only neutrons ( Z=0 and N=126) as shown in Fig. [8l As expected, the scaling 
function for the proton only results show a much larger dependence on the scattering angle 
due to the changing representation of the longitudinal and transverse response. The neutron 
only scaling function contains much less longitudinal contributions (due to the motion of 
the magnetic moment) and varies less with electron scattering angle. 

In Fig. [9|, we investigate scaling at larger values of the four momentum transfer as a func- 
tion of electron scattering angle 6 for quasielastic scattering on ^osp]^ when both Coulomb 
distortion and final state interaction are included. Based on this result and other cases 
we have examined, scaling is only weakly dependent on 6 at these higher kinematic values 
although there is still some deviation at very large electron scattering angles. Furthermore, 
as shown in Fig{TO] scaling for different nuclei still works well if all cases are calculated at 
the same backward scattering angle. Clearly as the final state interaction decreases scaling 
improves. In Fig. [HI we show the new scaling functions at high incident energy, E = 1.5 
GeV, for the scattering angle 6 = 30° in terms of y' as in Fig. [5l Clearly scaling is greatly 
improved even at this intermediate four momentum transfer region as long as we do not 
vary the electron scattering angles. 

In summary, we have investigated y-scaling of the inclusive (e, e') reaction from A=12, 
20, 56, 197, and 208 in the quasielastic region which is kinematically isolated from inelastic 
scattering. We use a realistic nuclear model describes the available data quite well and 
our calculations include electron Coulomb distortion. As shown in Figs. [6] and [HI the 
scaling function F{y) is not independent of electron scattering angle with particularly strong 
dependence at backward scattering angles for lower values of the four momentum transfer. 



However, we have demonstrated that if one restricts the variation of the electron scattering 
angles (thereby not changing the longitudinal transverse mix significantly) , scaling still holds 
at these lower kinematic values. Furthermore, as the energy increases, the dependence 
of scaling on electron scattering angles is greatly reduced. Wc also found that y-scaling 
breaks down at lower kinematic values due to the final state interaction. Again, as the 
energy increases (and the final state interaction weakens) the breakdown of scaling is greatly 
reduced. We also used our model to investigate the different binding energy effects and final 
state interaction effects in different target nuclei. We do find some improvement in scaling by 
modifying the definition of the scaling variable y to include binding effects. Finally, electron 
Coulomb distortion does disrupt scaling to some degree, but does not play a significant role. 

In conclusion, we confirm that the final state interaction is the primary cause of the scal- 
ing violation in the quasielastic region for lower energies, but we find that scaling improves 
dramatically at larger four momentum transfer. Our results show that at larger four mo- 
mentum transfer values, scaling holds across a range of nuclei quite well. Our model results, 
which agree very well with the available data, suggest that scaling can be used at larger 
four momentum transfer values to subtract out the quasielastic contributions to the (e, e') 
cross section from nuclei so that inelastic contributions which contribute incoherently to the 
measured cross section can be investigated. Furthermore, our results suggests that as the 
energies increase the scaling function obtained from different nuclei and at different electron 
scattering energies differ at most by only 10 to 20 %. 
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FIG. 1: The y-scaling functions for ^^C, ^''Ca, ^^Fe, ^^^Au, and ^''^Pb with incident electron energy 
750 MeV and scattering angle 45°. The calculations do not include the final state interaction of the 
outgoing nucleons nor electron Coulomb distortion at the four momentum squared of approximately 
0.2 (GeV/c)2 < Q2 < 0.3 (GeV/c)^. 




FIG. 2: The same as in Fig. [T] except electron Coulomb distortion is included. 
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FIG. 3: The same as in Fig. [T] except both electron Coulomb distortion and the final state 
interaction of the outgoing nucleons are included. 
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FIG. 4: The y-scaling functions for three different target nuclei from ^^C (sold and •), ^^Fe (dot 
and ■), and ^^"^Au (dash and A). The electron incident e nerg y is E = 2.02 GeV and the scattering 




angle is ^ = 15°. The experimental data are from SLAG 
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FIG. 5: The same kinematics as in Fig. [3] except plotted as a function of the new y'-scahng variable. 
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FIG. 6: The y-scaling functions for '^'^Ca at three different electron energies, E = 739 MeV and 
scattering angle 6 = 45.5° (sohd and •), E = 375 MeV and 9 = 90° (dot^and ■), and E = 367 
MeV and 9 = 140° (dash and A). The experimental data are from Bates 
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FIG. 7: The y-scaling functions for only the protons from ^'^^Pb at three different electron energies, 
E = 550 MeV and scattering angle 6 = 60°, E = 354 MeV and 9 = 90°, and E = 310 MeV and 
9 = 143° with Coulomb distortion and final state interaction included. 
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FIG. 8: The same as in Fig. [7| except for neutrons only. 
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FIG. 9: The y-scaling functions for ^''^Pb both with final state interaction and Coulomb distortion 
at five different electron kinematics, E = 1830 MeV and scattering angle 6 = 30° (solid and red), 
1000 MeV and 60° (dash and black), 740 MeV and 90° (dot and blue), 630 MeV and 120° (dash- 
dot and skyblue), and 580 MeV and 150° (two-dot and black). The range of the four momentum 
squared is approximately 0.5 (GeV/c)^ < < 0.8 (GeV/c)^. 
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FIG. 10: The 7/'-scaling functions for several nuclei with final state interaction and Coulomb 
distortion at the incident electron energy E = 580 MeV and scattering angle 9 = 150°. The range 
of the four momentum squared is approximately 0.5 (GeV/c)^ < < 0.8 (GeV/c)^. 
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FIG. 11: The new y'-scaling functions for the high electron energy E = 1.5 GeV and scattering 
angle 6 = 30° from several nuclei. 



